Astronomy & Astrophysics manuscript no. astroph 


© ESO 2009 


April 7, 2009 





o 
o 



O 



O 

o 



X 



A break in the gas and dust surface density of the disc around the 

T Tauri star IM Lup 

O. Panic^ M. R. Hogerheijde' , D. Wilner^, and C. Qi^ 



' Leiden Observatory, Leiden University, P.O.Box 9513, 2300 RA, Leiden, The Netherlands 

- Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 



Received ; Accepted 



ABSTRACT 



Aims. We study the distribution and physical properties of molecular gas in the disc around the T Tauri star IM Lup on scales close to 
200 AU. We investigate how well the gas and dust distributions compare and work towards a unified disc model that can explain both 
gas and dust emission. 

Methods. '^CO, '^CO, and C'**0 7=2-1 line emission, as well as the dust continuum at 1.3 mm, is observed at I'.'S resolution towards 
IM Lup using the Submillimeter Array. A detailed disc model based on the dust emission is tested against these observations with the 
aid of a molecular excitation and radiative transfer code. Apparent discrepancies between the gas and dust distribution are investigated 
by adopting simple modifications to the existing model. 

Results. The disc is seen at an inclination of 54° ±3° and is in Keplerian rotation around a 0.8-1.6 M© star. The outer disc radius 
traced by molecular gas emission is 900 AU, while the dust continuum emission and scattered light images limit the amount of dust 
present beyond 400 AU and are consistent with the existing model that assumes a 400 AU radius. Our observations require a drastic 
density decrease close to 400 AU with the vertical gas column density at 900 AU in the range of 5 x 10^"- 10"^ cm"^. We derive a 
gas-to-dust mass ratio of 100 or higher in disc regions beyond 400 AU. Within 400 AU from the star our observations are consistent 
with a gas-to-dust ratio of 100 but other values are not ruled out. 

Key words, planetary systems: protoplanetary discs - stars: individual: (IM Lup) - stars: pre-main-sequence 



1. Introduction 

Low-mass star formation is commonly accompanied by the for- 
mation of a circumstellar disc. The disc inherits the angular 
momentum and composition of the star's parent cloud, and is 
shaped by the accretion and other physical processes within 
the disc during the evolution that may result in a planetary 
system. Over the past two decades observations of circumstel- 
lar discs at millimetre wavelengths have been established as 
powerful probes of the bulk of the cold molecular gas and 
dust. Spatially resolved observations of the molecular gas with 
(sub) millimetre interferometers constrain the disc size and in- 
clination, the total amount of gas, its radial and vertical struc- 
ture, and the gas kinematics (e.g . , [G uilloteau & DutrevI [T998 

I4tllsel] 



llsella et al. 2007; 



Dartois et al. 2003; Oi et al.ll2004 

2007; Panic et al, 2008h . In parallel, continuum observations 
of the dust at near-infrared to millimetre wavelengths provide 
the disc spectral energy distribution (SED), that through mod- 
elling (e.g., [Ch iang & Goldreich 1997; Dullemond et al. 200jl 
iD'Alessio et al.ll200 5) yields the disc's density and temperature 
structure from the disc inner radius to a few hundred AU from 
the star. Studies of the gas through spatially resolved molecu- 
lar line observations using results from the SED modelling (e.g., 
[Raman et al. 2006; Panic et al. 2008) offer the means of address- 
ing the gas-to-dust ratio, differences between the radial and verti- 
cal distributions of the gas and the dust, and the thermal coupling 
between the gas and the dus t in the upper disc la yers exposed 
to the stellar radiation (e.g.. iJonkheid et aLll2004l) . Recent pa- 
pers have suggeste d different disc siz es for the dust and the gas 
(e.g., HD 163296. Ilseila et al.ll2007h . which may be explained 
by sensitivity effects in discs with tapered outer density pro- 



files (iH ughes et alj|2008l) . Here, we present the results of a com- 
bined study using spatially resolved molecular-line observations 
and SED modelling of the disc around the low-mass pre-main- 
sequence star IM Lup. 

Most pre-main-sequence stars with discs studied so far in 
detail are located in the nearby star-forming region of Taurus- 
Aurigae, accessible for the millimetre facilities in the north- 
ern hemisphere. Much less is known about discs in other star- 
forming regions such as Lupus, Ophiuchus or Chamaeleon. It 
is important to compare discs between different regions, to in- 
vestigate if and how different star-forming environments lead 
to differences in disc properties and the subsequent plane- 
tary systems that may result. IM Lup is a pr e-main-sequence 
star lo cated in the Lupus II cloud for which Wichmann et al.l 
(119981) report a distance of 190+27 pc using Hipparcos par- 
allaxes. From its MO spec tral type and estimated bolometric 
luminosity of 1.3+0.3 L©, Hughes et al.l ( fl994 ) derive a mass 
of 0. 4 Mq and an age of 0.6 Myr using evolutionary tracks 
from ISwenson eta U (1994) , or 0.3 M p and 0. 1 Myr adopt- 
i ng th e tracks of [P'Antona & Mazzit elli (1994). In Pinte et al| 
(2008 ), a much high er value of 1 M© is derived using tracks of 
Baraffeetal.l(ll998h . 

IM Lup is surrounded by a di sc detected in sca ttered light 
with the Hubble Space Telescope (iPinte et al.l 12008') and in the 
CO 7=3-2 line w i th the James Clerk Maxwel l Telescope by 
Ivan Kempen et al.l (l2007l) . iLommen et al.l (|2007|) conclude that 
grain growth up to millimetre sizes has occur ed from con t inuum 
measurements at 1.4 and 3.3 mm. Recently. i Pinte et al] (|2008|) 
present a detailed model for the disc around IM Lup based on 
the full SED, scattered light images at multiple wavelengths 
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from the Hubble Space Telescope, near- and mid-infrared spec- 
troscopy from the Spitzer Space Telescope, and continuum 
imaging at 1 .3 mm with the Submillimeter Array. They conclude 
that the disc is relatively massive, M w 0.1 with an uncer- 
tainty by a factor of a few, has an outer dust radius not greater 
than «400 AU set by the dark lane and lower reflection lobe seen 
in the scattered light images, and has a surface density pro- 
portional to ' constrained by the 1.3 mm data. Furthermore, 
they find evidence for vertical settling of larger grains by com- 
paring the short- wavelength scattering properti es to the grain- 
size co nstraints obtained at 1.4 and 3.3 mm bv iLommen et al.l 
(I2007h . 

In this work, we present (Sect. |2]i spatially resolved obser- 
vations of the disc around IM Lup in '^CO, '^CO and C"*0 
7=2-1 line emission, together with 1.3 mm dust continuum 
data, obtained with the Submillimeter Array{3 (SMA). Our re- 
sults (Sect. 13) show that the gas disc exten ds to a radius o f 
900 AU, more than twice the size inferred by P inte et al](l2008l). 
A deta iled comparison (Sect. 14. Il l to the model of Pinte et al] 
(1200 8') suggests a significant break in the surface density of both 
the gas and the dust around 400 AU, and we discuss possible 
explanations. We summarise our main conclusions in Sect. |5] 



2. Observations 

IM Lup was observed with the SMA on 2006 May 21 in a 
8.6-hour track, with a 4.3 hours on-source integration time. 
The coordinates of the phase centre are RA= 15*'56™09U7 and 
Dec= -37°56'06'.'40 (J2000). Eight antennas were included in 
an extended configuration providing a range of projected base- 
lines of 7 to 140 meters. The primary beam half-power width is 
55'.'0. The SMA receivers operate in a double-sideband (DSB) 
mode with an intermediate frequency band of 4-6 GHz which 
is sent over fiber optic transmission lines to 24 "overlapping" 
digital correlator chunks covering a 2 GHz spectral window in 
each sideband. The DSB system temperatures ranged from 90 to 
150 K. The correlator was configured to include the '^CO 7=2- 
1 line (230.5380000 GHz) in the upper sideband and the '^^CO 
2-1 (220.3986765 GHz) and C'-^O 2-1 line (219.5603568 GHz) 
in the lower sideband. Each of the three lines was recorded in a 
spectral band consisting of 512 channels with 0.2 MHz spac- 
ing (~0.26 km s"'). Simultaneously to the line observations, 
the 1.3 mm dust continuum was recorded over a bandwidth of 
1.8 GHz. 

The data were calibrated and edited with the IDL-based 
MIR software package (0). The bandpass response was deter- 
mined from Jupiter, Callisto and 3C273. After the passband 
calibration, broadband continuum in each sideband was gener- 
ated by averaging the central 82 MHz in all line-free chunks. 
Complex gain calibration was performed using the quasar 
J1626-298. The absolute flux scale was set using observations 
of Callisto. Subsequent data reduction and i mage analysis wa s 
carried out with the Miriad software package ( Sault et akll 19951) . 
The visibilities were Fourier transformed with natural weight- 
ing, resulting in a synthesized beam of r.'8 x 1'.'2 at a position 
angle of 0.2°. 1 Jy/beam corresponds to 15.9 K. The r.m.s noise 
level is 125, 94 and 102 mJy beam ' per channel respectively for 
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' The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 

^ http://www.cfa.harvard.edu/~cqi/mircook.html 



Fig. 1. Dust continuum image at 1.3 mm. The contours are at (2, 
4, 8, 16)x3.67 mJy beam"' (2, 4, 8, 16 sigma) levels. The filled 
ellipse in the lower left corner shows the synthesized beam. 

the '-CO, '^CO and C'^O spectral line data and 4 mJy beam"' 
for the continuum data. 



3. Results 

3.1. Dust Continuum 

Figure[T]shows the 1.3 mm continuum emission observed toward 
IM Lup, previously reported in Pi nte et al. (2008). The emission 
is centered on RA=15''56"09^20, Dec= -37°56'06'.'5 (J2000), 
offset by -i-0'/4 in right ascension and by -O'.'l in declination 
from the pointing center We adopt the peak of the continuum 
emission as the position of IM Lup. The peak intensity of the 
continuum emission is 104 + 4 mJy beam"' and the total flux 
176+4 mJy. The emission intensity is fit to the precision of one 
sigma by an elliptical Gaussian, yielding a source FWHM size 
of 1'.'52+0'.'15 X l'.'06+a.'15 and aposition angle of -35°5 +4?0 
deconvolved with the synthesized beam. This position angle, and 
the inclination in the range of 33°-53° suggested by the de- 
convolved a sgect ratio, agree well with the values obtained by 
iPinte et all (|2008') of, respectively, -37°+5° and 45°-53° from 
scattered light imaging. 

A fit to the 1 .3 mm visibilities done in lPinte et al.l (l2008l) pro- 
vides a rough disc mass estimate of 0. 1 M©, with an uncertainty 
of a factor of few, dominated by the adopted dust emissivity and 
gas-to-dust mass ratio in the model. 



3.2. Molecular Lines 

Emission of '^CO and '^CO 7=2-1 was detected toward IM Lup, 
and an upper limit on C'^O 2-1 obtained. Figure |2] shows the 
zero moment (integrated emission, contours) and first moment 
(velocity centroid, colour scale) of the '^CO and '^CO emission 
from IM Lup. Significantly detected '^CO emission extends to 
5" from the star (roughly 900 AU for IM Lup). This is more 
than double the size inferred from the dust continuum image, and 
Sect.|4]discusses if this is due to different sensitivity in these two 
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Fig. 2. (a),(b): First moment maps in the '^CO and '^CO 7=2-1 
lines, from 1.9 kms"' to 6.9 kms"' observed towards IM Lup. 
These maps are created using the Miriad task 'moment' with 
clip levels of 0.5 and 0.3 Jy respectively. The integrated emis- 
sion of '^CO 7=2-1 is shown in contours of 1, 2, 3, ...x500 mJy, 
and that of '^CO 7=2-1 with 1, 2, 3, ...xl60 mJy contours, (c), 
(d): First m oment and integrated emission maps calculated using 
iPinte et aTl model and same clip level, velocity range and con- 
tour levels as in (a) and (b). (e), (f): First moment and integrated 
emission maps calculated using extended disc model (described 
in Sect. |4.2l l with model parameters Sigma4oo =2x10^' cm"^ 
and p =1. The clip level, velocity range and contour levels are 
as in (a), (b), (c) and (d). 



tracers or if the gas disc indeed extends further than the dust disc. 
The aspect ratio (5/3), suggesting an inclination of 53° + 4°, and 
orientation PA=-36+5° of the CO disc agree with the continuum 
image (Sect. [3Tt and scattered light imaging results (Pinte et al] 
120081) . 

The first moment images of Fig. |2] show velocity patterns 
indicative of a rotating disc inclined with respect to the line 
of sight. This is also seen in Fig. [3] which presents the '^CO, 
'^CO, and C"^0 spectra averaged over 8" x 8" boxes around IM 
Lup. The '^CO and '-^CO lines are double-peaked and centered 
on Vlsr =4.4+0.3 km s Figures |4] and |5] show the channel 
maps of the '^CO and '^CO emission, respectively, revealing the 
same velocity pattern also seen from the first-moment maps and 
the spectra. The Keplerian nature of the velocity pattern is most 
clearly revealed by Fig.|6] which shows the position-velocity di- 
agram of the '^CO emission along the major axis of the disc. In 
Section|4] we derive a stellar mass of 1.2 M©, and, as an illustra- 
tion, the rotation curves for stellar masses of 0.8, 1 .2, and 1 .6 M© 
are plotted in Fig. |6] 
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Fig. 3. '^CO, '^CO, and C'**0 7=2-1 line spectra summed over 
8" X 8" regions centered on the location of IM Lup. The '^CO 
and '^CO spectra are shifted upward by 2 and 5 Jy, respectively. 
The dashed red line shows the line centre at vl5«=4.4 km s '. 
The grey zone indicates the range from 4 to 6 km s"' where the 
'^CO line is significantly affected by the foreground absorption; 
the corresponding part of the '^CO spectrum is plotted with a 
dotted line. 



Us ing single-dish '^CO 3-2 observations. 'van Kempe n et al.l 
(l2007h first identified molecular gas directly associated with 
IM Lup, but they also conclude that the vis R-range of 4 to 
6 km s"' is dominated by gas distributed over a larger (> 
30") scale. In our '^CO 2-1 data this same visR-range is 
also likely affected: wh ere the single-dish '"CO 3-2 spectrum 
from Ivan Kempen et al.l shows excess emission over visr =4- 
6 km s~', the red peak of our '^CO 2-1 spectrum, which lies 
in this VLSR-iange, is weaker than the blue peak at -1-3.5 km s 
We susp ect that absorptio n by the same foreground layer identi- 
fied by v an Kempen et a l.'is resonsible for this decrement, while 
its emission is filtered out by the interferometer The '^CO 2-1 
spectrum is symmetric, suggesting that the foreground layer is 
optically thin in this line. 

The spatial extent of the line emission is further explored in 
Fig. |7] which plots the '^CO and '^CO 7=2-1 vector -averaged 
line fluxes against projected baseline length. The '^CO flux is 
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Fig. 4. The black contours show the observed '^CO J -2-1 emission in the velocity range from 2.46 to 6.42 km s"'. Alongside the 
observations, the panels with grey contours show the calculated emission from the extended disc model described in Sect. 14.21 with 
parameters E400 —2x lO^^cm"^ and p —I. Labels indicate the velocity of each channel. The lower left corner of bottom-left panel 
shows the size and position angle of the synthesized beam. The contour levels are -1, 1, 2, 3, 4x400 mJy beam"' (~3cr) in all 
panels. 
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Fig. 6. Position-velocity diagram of the i^co 2-1 hne erms- 
sion along the disc's major axis. Contour levels are (1, 2, 3, 
...)x400 mJy (~3cr). For comparison, the thick solid line cor- 
responds to Keplerian rotation around a 1.2 M© star; dashed 
and dotted lines correspond to the stellar masses of 0.8 and of 
1 .6 Mq, respectively. 



integrated from 2.5 to 4.0 kms ' to avoid the range where fore- 



ground absorption affects the line. The '^CO flux does not suffer 
from absorption and is integrated over its full extent from 2.5 to 
6.9 kms"' . Comparing the curves of Fig.|7]to those of the contin- 
uum flux versus baseline lengths (Fig.[8]l it is clear that the line 
flux is much more dominated by short spacings (<40 k/l). This 
profile may indicate the presence of a larger structural compo- 
nent (outer disc or envelop e), combined with the disc emission 
(See |j0rgensen et aLll2005L Fig. 2). We explore disc structure 
beyond 400 AU in the following section. 



4. Discussion 

The results of the previous section show that IM Lup is sur- 
rounded by a gaseous disc in (roughly) Keplerian rotation. The 
gas disc has a radius of 900 AU, and its surface density may 
have a break around a radius of 400 AU. In contrast, the size 
of the dust disc is constrained to a r adius of 40 AU b y our 
continuum data and the modelling of iPinte et alj ( l2008l) . This 
section explores if the gas and dust have the same spatial dis- 
tribution (in which case different sensitivity levels need to ex- 
plain the apparent difference in size) or if the gas and dust are 
differentl y distributed radia lly. First we investigate whether the 
model of Pinte et al. I (I2008h can explain the molecular line ob- 
servations (Sect. |4~TT l. After we conclude that this is not the case, 
we construct new models for the gas disc (Sect. |4.2| i describ- 
ing their best-fit parameters, and compare them to the dust disc 
(Sect. 1431). 
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Fig. 5. Channel maps of the observed '^CO 7=2-1 emission at the velocities where the line is detected are shown in black contours. 
For comparison, the line emission calculated from our extended disc model described in Sect. 14.21 is shown in grey contours. The 
model parameters are S400 = 2 x 10^'cm"^ and p -I. Labels indicate the velocity of each channel. The lower left corner of bottom- 
left panel shows the size and position angle of the synthesized beam. The contour levels are -1, 1, 2, 3, 4x200 mJy beam"' (~2cr) 
in all panels. 




Fig. 7. The left and the right panel show vector-averaged CO and CO line flux (black symbols), respecti vely. Dashed black lines 
represent the zero-signal expectation value of our line visibility data. The calculated visibilities based on ' Pinte et aU model (full 
red line) and our extended disc model described in Sect. 14.21 (dotted red line) are shown for comparison. Our model parameters are 
£400 = 2 X 10^' cm"^ and p =1. The '^CO flux is integrated over the 0.8-4.0 kms"' range and '''CO over 2.5-6.9 kms"', covering 
the full line width. 



4.1. Molecular-line emission from the dust-disc model 

Recently, iPinte et akl ( l2008 l) present a detailed model for 
the disc around IM Lup based on the full SED, scattered 
light images at multiple wavelengths from the Hubble Space 
Telescope, near- and mid-infrared spectroscopy from the Spitzer 



Space Telescope, and continuum imaging at 1.3 mm with the 
Submillimeter Array. 

Based on the two-dimensional density and temperature 
structure of the Pinte et al. model, with M = 0.1 Mq, R„„i = 
400 AU, and / = 50°, we calculate the resulting line intensity of 
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Fig. 8. Vector-averaged continuum flux as a function of projected 
baseline length (black symbols). Error bars show the variance 
within each annular bin. The dashed histogram shows the zero- 
signal expectation v alue. The fu ll red line shows the continuum 
flux calculated from lPinte et al.l model. The dotted red line cor- 
responds to the extended disc model with ^400 - 2x10^' cm"^ 
and p -I (see Sect. 14. 2l for model description). 

the '^CO and '^CO 7=2-1 lines. To generate the input model 
for the molecular excitation calculations, we adopt a gas-to- 
dust mass ratio of 100 and molecular abundances typical for the 
dense interstellar medium ( Frerking et al.|[T982HWilson & RoodI 
11994 : a '^CO abundance with respect to H2 of 10"'* and a 
'^CO/'^CO abundance ratio of 77. No freeze-out or photodis- 
sociation of CO is included. The velocity of the material in the 
disc is described by Keplerian rotation around a 1.0 star 
plus a Gaussian microturbulent velocity field with a FWHM of 
0. 16 km s ' ; the exact value of the latter parameter has little ef- 
fect on the results. Using the molecular excitation and radi ative 
transfer code RATRAN dHogerheiide & van der 'fag '2000') and 
CO-H2 collision rates fro m the Leiden A tomic and Molecular 
Database (lambdaB Sch oier et al.l l2005) we calculate the sky 
brightness distribution of the disc in the '^CO and '^CO 7=2-1 
lines for its distance of 190 pc. From the resulting image cube, 
synthetic visibilities corresponding to the actual (m, v) positions 
of our SMA data were produced using the MIRJAD package 
dSault et al.lll995h . Subsequent Fourier transforming, cleaning, 
and image restoration was performed with the same software. 

Figure |2] compares the zeroth-moment (integrated intensity; 
contours) and first-moment (velocity-integrated intensity; colour 
scale) maps of th e resulting synthetic images to the observa- 
tions. Clearly, the iPinte et al.l model produces '^CO and '^CO 
2-1 emission with spatial extents and intensities too small by a 
factor close to two. In Fig.|7]it is clear that the Pinte et al. model 
fails to reproduce the '^CO and '^CO line fluxes at short pro- 
jected baseline lengths, but is consistent with the observations 
longward of 40 kA th at correspo nd to spatial scales <500 AU. 
Our comparison with Pinte et al. model thus suggests that the 
gas extends much further than 400 AU from the star. 

The observed 1.3 mm continuum emission tra ces the ex- 
tent of larger dust particles (up to millimetre in size). |Pinte et al.l 
(I2OO8I) show that their 400 AU model reproduces these observa- 
tions well. In Sect. 14. 31 we explore to what level larger particles 
can be present outside 400 AU. 

' http://www.strw.leidenuniv.nl/~moldata/ 



4.2. Extending the gas disc beyond 400 AU 

As mentioned in Sect. 13.21 the CO Une flux as function of pro- 
jected baseline length suggests a possible break in the emis- 
sion around40 kA (Fig. |7]l. Results of Sect. 14.11 show that 
the iPinte et akl model, while providing a good description of 
line fluxes at small spatial scales (baselines> 4QkA), requires 
a more extended component to match the observed line fluxes 
(baselines< 40k/l). In this section we extend the iPinte et aTl 
model by simple radial power laws for the gas surface density 
and temperature, and place limits on the gas column densities in 
the region between 400 and 900 AU. 

Table [T] lists the model parameters. For radii smaller than 
400 AU, the radi al and verti cal density distribution of the ma- 
terial follows the iPinte et aLl model. As in Sect. 14.11 we adopt 
'standard' values of gas-to-dust mass ratio and molecular abun- 
dances, and a Gaussian microturbulent velocity field with equiv- 
alent line width of 0.16 km s '. Unlike the calculations of 
Sect. 14.11 we add as free parameters the stellar mass and 
the gas kinetic temperature. For the latter, we fo llow the two- 
dimensional structure prescribed by iPinte et all but scale the 
temperatures upward by a factor / with 1 < / < 2. This cor- 
responds to a decoupling of the gas and dust temperatures, as 
may be expected at the significant height above the midplane 
where t he '^CO and '^CO lines originate (see, e.g., lOi et al.l 
I2OO6; Jo nkheid et"aLl |2004 . Because the highly red- and blue- 
shifted line emission (line wings) comes from regions closer to 
the star than 400 AU and is optically thick, factor f is determined 
by the observed fluxes in the line wings. The molecular excita- 
tion and synthetic line data are produced in the same way as 
described in Sect. 14. II 

Outside 400 AU we extend the disc to 900 AU, as sug- 
gested by the observed '^CO image of Fig. |2l by simple 
radial power laws for the surface density and temperature, 
2=2400 (R/400 AU)-'' and T = ^400 (R/400 AUy^. At 400 AU, 
the surface density is ^400 and the temperature is 7400; the pa- 
rameter p is assumed to be > 0. To limit the number of free 
parameters, we set 7400 = 30 K and q = 0.5; we assume that the 
disc is vertically isothermal and that the '^CO abundance is lO"'*, 
constant throughout the disc. At R >400 AU, the disc thick- 
ness is set to Zmax=100 AU and the density p{R, z)-'^(R) / Zmax 
is vertically constant. For our free parameters £400 and p, we 
assume that £400 <0.9 g cm"^ (vertical gas column density of 
2x10^^' cm"^), the value at the outer radius of the iPinte et al.l 
model. We have run a number of d isc models, with the inner 
400 AU described by the lPinte et alj model (with the gas kinetic 
temperature scaled as described in the previous paragraph) and 
the region from 400 to 900 AU described here by the disc ex- 
tension. Figure [To] shows the surface density in the models that 
we have t ested: within 400 AU it is the surface density as in 
iPinte et al . (blue line) and between 400 and 900 AU different 
combinations of E400 and p (black lines). The models are tested 
against the observed '^CO and '^CO Mv-data, channel maps, 
spectra, and position-velocity plots. The comparison of mod- 
elled emission with wv-data for the line wings, vlsr <3.0 km s ' 
for '^CO and vlsr <3.5 km s"'-i-vlsr >5.5 km s"' for '•'CO is 
also examined. 

Figure [10] shows the models that overproduce the observed 
emission with dashed black lines and those that underproduce it 
with dotted black lines. The full black lines correspond to the 
models that reproduce well our '^CO and '^CO data. The gen- 
eral area (beyond 400 AU) allowed by the models is shaded in 
Fig.[TO]for guidance. It can be seen that the '^CO and '^CO ob- 
servations constrain the column density of '^CO atR- 900 AU 
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Table 1. Model parameters 



Parameter 


P ^AOn ATT 


400< R <800 AU 


Surface density 


E oc R-\ see Pinte et al. (2008) 


S400 (/?/400 AU)-" with E400 <0.9 g cm" 


^ P>0 


Gas-to-dust mass ratio 


1 nn 


100 




Gas temperature structure 


/Td^s, withl </<2 


7400 (i?/400 AU)-", 7400 = 30 K, 9 = 


-- 0.5 


Vertical structure 


see Pinte et al. (2008) 


T{z) = constant, p(z) = constant, Zma\ = 


100 AU 


[CO]/[H2] 


lO'' 


10-* 




['-C0]/['3C0] 


77 


77 






1.2 Mo 


1.2 Mo 




Inclination 


50° 


50° 




FWHM microtubulence 


0.16 km S-' 


0.16 km S-' 





toA^co = (0.05 -1.0)x 10'^ cm-^, where the lower bound follows 
from the requirement that the '^CO emission is sufficently ex- 
tended and the upper bound from the requirement that the '^CO 
and '-^CO peak intensity, and the extent of the '^CO emission are 
not overestimated. The corresponding surface density at 900 AU 
is Sgoo = (0.2 - 4.0) X 10-^ g cm-^, i.e., a vertical gas column 
density (0.05-1 .0) x 10^^ cm-^. Our data do not constrain the pa- 
rameters 2400 and p, that determine how the surface den sity de- 
creases from its value at the outer edge of the lPinte et al.l model, 
to its value at 900 AU. This is either a marked change from the 
power-law slope of p = 1 found inside 400 AU to p = 5 beyond 
400 AU, or a discontinuous drop by a factor ~1 0-1 00 in surface 
density at 400 AU. 

Figure |7] compares observations to synthetic '^CO and '-^CO 
line visibilities for our model with S400 = 2 x 10^'cm-^ and 
p - I, plotted with dotted red lines. There is a good match be- 
tween the model and the data for both transitions. In particular, 
the model reproduces well the change in the slope of visibili- 
ties, mentioned in Sect. 13.21 The match between the model (red 
lines) and observations is also seen in the line spectra. Fig. |9l 
Figures |4] and |5] show, respectively, the '^CO and '^CO chan- 
nel maps (black contours) compared to our extended disc model 
(grey contours). It can therefore be seen that our model provides 
a good description not only of the line intensity at each channel 
(spectra), but also a very close match in the spatial extent of the 
emission in each spectral channel. 

A good fit to the wings of the '^CO and '^CO spectra (Fig.|9l 
red lines) and the spatial distribution of the respective line fluxes 
at highly blue- and red-shifted velocities (lower panels, Figs|4] 
and|5]) is found for temperature scalings / of 1.7 for '^CO and 
1 .4 for '^'CO. These values of / suggest that the gas is somewhat 
warmer than the dust at the heights above the disc where the 
'^CO and '^CO emission originates, and more so at the larger 
height probed by the ^^CO line compared to the '^CO line. 

At the adopted disc inclination of 50°, the line peak separa- 
tion provides a reliable contraint on the stellar mass. We find a 
best-fit of M* = 1 .2 ± 0.4 Mq, where the uncertainty is domi- 
nated by our limited spectral resolutio n. This value is con sistent 
with the rough estimate of 1 from Pinte et al.' (2008), but a 
few times higher than derived by Hughes et al. ( 1994). 

We conclude that the surface density traced through ^^CO 
and '^CO has a discontinuity around R =400 AU either in 
^co{R) or in its derivative dYcoldR, or both. This may, or may 
not be an indication of an overall discontinuity of the gas sur- 
face density. A break in the temperature T{R) cannot explain 
the observations, since our model akeady adopts a low temper- 
ature at the margin of '^CO freeze-out in the outer regions. An 
alternative explanation for the observations is a radical drop in 
the abundance of CO (with respect to H2 and H) or its radial 
derivative. Freeze out onto dust grains or photodissociation can 
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Fig. 9. '^CO and '^CO 7=2-1 Hne spectra avera ged over a S'.'O x 
S'.'O region centered on IM Lup. The '^CO spectrum is shifted 
up by 3 Jy for clarity. The red lines show the emission predicted 
the extended disc model described in Sect. 14. 21 with E400 = 2 x 
10^'cm-^ and p - I. 



significantly reduce the gas-phase abundance of CO. In the next 
section we explore the limits that the dust emission can give us 
on the gas content outside 400 AU, and compare them to the 
'^CO results. 
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Fig. 10. Gas surface density in our disc models is plotted as a 
function of radius. Within 400 AU, it is identical to the Pinte 
et al. model shown with the full blue line. Outside 400 AU, 
we explore different power-law distributions, each plotted in 
black and marked with the corresponding slope p. The models 
which overestimate the observed '^CO emission are plotted with 
dashed lines, while those that underpredict it are shown in dot- 
ted lines. The full black lines represent the models that fit well 
the '^CO 7=2-1 emission, and define the shaded region which 
shows our constraint on Sco/[CO] in the outer disc, where a CO 
abundance of [CO] = 10"^ is used. The upper limit on Sdust^/i^ 
placed by scattered light images is shown with a dark grey line, 
with a gas-to-dust mass ratio ^/ii=100. For comparison, the red 
lines correspon d to disc mod e ls wit h an exponential drop off 
as described in iHughes et al.l (I2008h . The model shown with 
the long-dashed red line has parameters y =0.3, ci=340 AU 
and C2 - 3.1 x 10^"^ cm"- and fits the gas constraints. The 
model with the dash-dotted red line y =0.6, ci=340 AU and 
C2 = 1.8 X lO^-' cm"^ fits the scattered light constraints. No sin- 
gle model with a tapered outer edge can fit both these constraints 
and the constraints within 400 AU simultaneously. 

4.3. Comparing gas and dust at radii beyond 400 AU 

The previous section concluded that both the gas and the dust 
out to 400 A U in the di sc around IM Lup is well described by 
the model of Pinte et all with the exception of gas temperatures 
that exceed the dust temperature at some height above the disc 
midplane. It also found that the gas disc needs to be extended to 
an outer radius of 900 AU, albeit with a significant decrease in 
the surface density of CO, "Leo, or in its first derivative, c/2co / dR 
clo se to 400 AU 

iPinte et akl ( |2008|) show that some dust is present outside 
400 AU as well, visible as an extended nebulosity in their 0.8 fim 
scattered light images. At the same time, the visible lower scat- 
tering disc surface places a stringent limit on the surface den- 
sity Sdust of small dust particles outside 400 AU. Requiring 
the optical depth t = 'LdustK < 1 and adopting an emissiv- 
ity per gram of d ust of /c=(8-10)xI0^ cm^g~' at 0.8 fim (See 
first row of Tab. l. lOssenkopf & Henning|[T994h . we find S^„j, < 



(1.0 - 1.3) X lO"'* g cm"^. If we adopt the gas-to-dust mass ratio 
of 100, this corresponds to A^h, < (2.5 - 3 - 1) x 1 0^' cm"^. Our 
limit differs from that given in iPinte et akl ( |2008|) (0.2 g cm"^) 
because we use dust opacities representative of small dust, while 
they assume considerable grain growth in disc midplane and thus 
use much lower dust opacities at 0.8 fim. The limit on surface 
density we derive is two orders of magnitude low er than the 
column density at the outer radius of 400 AU of the i Pinte et aTl 
model. This indicates that either the dust surface density drops 
sharply at 400 AU, or that efficient grain growth beyond 400 AU 
has caused a significant decrease in dust near-IR opacity. As 
can be seen in Fig. [TO] the upper limit on surface density of 
(2.5 - 3.1) X 10^' cm"- is consistent with the gas surface density 
range inferred in Sect. l4.2l from our CO data, using the canonical 
CO/H2 abundance of lO l 

While 0.8 //m imaging traces the small dust, our observations 
of 1.3 mm dust continuum emission, on the other hand, trace 
the milhm etre-sized dust particles. In Fig. |8]we can see that the 
iPinte et al.i model (full red line), with the radius of 400 AU, com- 
pares well to the observed continuum flux at all projected base- 
line lengths. On the other hand, the comparison of the 1.3 mm 
visibilities to our extended disc model with 24oo=2xlO^' cm"^ 
and 7? =1 shows that the model overestimates emission at short 
Mv-distances (large spatial scales). A constant dus t emissiyity of 
2.0 cm g" (emissivity of mm-sized grains, as in lDrainell2006h 
was used throughout the disc in the calculation of 1 .3 mm fluxes. 
Our model indicates that any dust present in the outer disc re- 
gions must be poor in mm-sized grains, i.e., have low millimetre 
wavelength opacities, while dust within 400 AU has likely un- 
dergone grain growth as found bv IPinte et all ( l2008l) . This fur- 
ther supports our choice of k at 0.8 fim when estimating the up- 
per limit on dust column (see above). Therefore, a viable model 
for the disc of IM Lup c onsists of a n 'inner' disc extending to 
400 AU as described in iPinte et aTl augmented with an 'outer 
disc' extending from 400 to 900 AU with a significantly reduced 
surface density (with negligible mass) but standard gas-to-dust 
mass ratio and CO-to-H2 ratios. The S EP of this new model 
should not differ significantly to that of iPinte et al.l model, and 
is therefore expected to provide a good match to the observed 
SED oflM Lup. 

iHughes et al.l (l2008l) find that the apparent difference be- 
tween the extent of submillimetre dust and gas emission in sev- 
eral circumstellar discs can be explained by an exponential drop 
off of surface density which naturally arises at the outer edge of 
a viscous disc. In Fig. [10] we show how, with a careful choice 
of parameters (y =0.3, ci=340 A U and C2 = 3.1 x lO^"* cm"^), 
the model of Hug hes et al.l (|2008) (red long-dashed line) can re- 
produce the surface density distribution of the models which de- 
scribe well the '^CO 2-1 line emission. This model, and the one 
discussed below, are only examples. A proper modelling of IM 
Lup in the context of viscous disc models would require a re- 
vision of the entire disc structure both in terms of temperature 
and density, which is outside of the scope of the current work. 
We notice that the iHughes et al.l models cannot simultaneously 
comply with the gas and dust constraints in the outer disc and the 
surface density derived by Pinte et al. (2008) in the inner disc. 
This is illustrated by the lHughes et al.l (12008 ) model with param- 
eters y =0.6, ci=340 AU and C2 = 1.8 x 10^"* cm"^, shown with 
the red dash-dotted line in Fig. [TO] The surface density of this 
model outside 400 AU is in agreement with observational con- 
straints from gas and dust, but it is ro ughly two orders o f mag- 
nitude lower than suface density from lPinte et al.l (l2008l) within 
400 AU. 
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In the standard theory of viscous discs (See lPrin^lelll98ll) . 
irrespective of the initial density distribution, a radially smooth 
surface density distribution with a tapered outer edge is rapidly 
reached. If there is a significant change in the nature of the vis- 
cosity inside and outside of 400 AU, discontinuities in the equi- 
librium surface density may follow. Such changes could, for ex- 
ample, result from differences in the ionization structure of the 
disc or from a drop of the surface density below some critical 
level. Here we explore some scenarios that could explain this: 

A young disc. An extreme example of such a configuration is 
a disc where the inner 400 AU follows the standard picture of a 
viscous accretion disc, but where the region outside 400 AU has 
not (yet) interacted viscously with the inner disc. This outer re- 
gion may be the remnant of the flattened, rotating prestellar core 
that has not yet made it onto the viscous inner disc. This config- 
uration, reminiscen t of the material around the object L1489 IRS 
dBrinch et al.ll200'f ). suggests that IM Lup would only recently 
have emerged from the embedded phase. L1489 IRS showed 
clear inward motion in its rotating envelope. Our observations 
limit any radial motions in the gas between 400 and 900 AU to 
< 0.2 km s or 20% of the Keplerian orbital velocities at these 
radii. Furthermore, for the 900 AU structure to survive for the 
lifetime of IM Lup of 0.1-0.6 Myr, inward motions cannot ex- 
ceed ~I0~^ km s"'. Any mass inflow is therefore small, and the 
material between 400 and 900 AU is likely on Keplerian orbits. 

A companion body. Another explanation for the break in the 
disc density structure around 400 AU would be the presence 
of a companion near this radius. A companion of ~ 1 Mjup at 
400 AU could open a gap in the disc and affect the viscous disc 
spreading. No co mpanions at this separation are visible in the 
HST images of (Pinte et al.li200 8) or in A'-band direct imaging 
dGhez et al.lll997h . Whether these observations exclude this sce- 
nario is unclear: it requires modelling of the orbital evolution of 
a companion in a viscously spreading disc and calculation of the 
observational mass limits at the age of IM Lup. This is beyond 
the scope of this paper. 

Gas to dust ratio. While our model is consistent with stan- 
dard gas-to-dust and CO-to-H2 ratios beyond 400 AU, this is 
not the only solution. Instead of adopting these standard ra- 
tios, which requires explaining the drop in S or dI,/dR around 
400 AU, we can hypothesize that the gas (H2 or H) surface den- 
sity is continuous out to 900 AU and that both the CO-to-H2 and 
dust-to-gas ratios show a break around 400 AU. This scenario 
requires a drop between 400 and 900 AU of the CO abundance 
by a factor between 10 and 200, and of the dust-to-gas mass ratio 
by a factor >90. These drops can be sudden, with a discontinu- 
ity at 400 AU, or more gradual, with a rapid decline of the two 
ratios from 400 to 900 AU. Since a low amount of dust emission 
outside 400 AU is observed both at wavelengths of ~ 1 mm (our 
data) and ~ 1 /im dPinte et al.ll2008h . the overall dust-to-gas ratio 
is likely affected, and not just the individual populations of small 
and large grains. 

Dust radial drift and photoevaporation. If a large fraction 
of the dust is removed from the disc regions outside 400 AU, 
the increased penetration of ultraviolet radiation could explain 
the drop in '^CO surface density through increased photodisso- 
ciation (van ZadelhoflF et al. 2003). Radi a l drift o f dust particles 
due to the gas drag force d Whipple! 1 9721: 1 Weidenschiliin^l 9771) 
is a possible scenario in circumstellar discs. The difference in 
velocity between the dust, in Keplerian rotation, and gas, sub- 
Keplerian because of the radial pressure gradient, can cause dust 
particles to lose angular momentum and drift inward. The opti- 
mal drift particle size depends on the gas density, Keplerian rota- 
tion frequency and hydrostatic sound speed. Most dust evolution 



models focus on the inner 100 AU of discs, relevant to planet for- 
mation. In these regions, the grains from 100 to about 0.1 m 
efficiently migrate inwards on a timescale shorter than 2 Myr 
However, the optimal grain size for inward drift decreases with 
the gas density. Our modelling of the disc region from 400 AU 
to 900 AU, predicts surface densities of ~ lO^^* g cm"^, low 
enough even for sub-micron-sized particles to drift inward (to 
< 400 AU). For the estimated age of IM Lup of 0.1-0.6 Myr, 
all particles larger than 0.1-0.02 /im will have migrated inward. 
This process leaves the outer disc unshielded by dust against UV 
radiation. Infrared emission of PAHs may be used to tr ace the 
disc surface in this scenario. However. iGeers et al.l d2007) do not 
detect PAH emission at 3.3 pm in their VLT-ISAAC L-band ob- 
servations of IM Lup. This may indicate that either there are not 
enough PAHs in the disc or that they are not exposed to a sig- 
nificant level of UV flux. The latter possibility allows the outer 
disc to remain molecular. Otherwise, the outer disc is exposed 
to photodissociating radiation, destroying much of the CO and 
likely also a significant fraction of the H2 given the limit on the 
dust surface density of 10^' cm"^ corresponding to Ay ^ 1'™^. In 
this scenario, the outer disc between 400 and 900 AU would be 
largely atomic and possibly detectable through 21 cm observa- 
tions of H I, or line observations of C I at 609 and 370 pm or C II 
at 158 pm. If photoevaporation is efficient in this region it may 
remove the (atomic) gas and reduce the gas surface density fur- 
ther Therefore, a combined effect of efficient drift, photodisso- 
ciation and photoevaporation in the outermost disc regions may 
be a reason for the low gas and dust density observed. The effi- 
ciency of these processes decreases with density and perhaps the 
density at 400 AU is high enough so that material is no longer 
efficiently removed from the disc. Only the detailed simultane- 
ous modelling of drift, photodissociation and photoevaporation 
could test this scenario. 

5. Conclusions 

We probe the kinematics and the distribution of the gas and dust 
in the disc around IM Lup through molecular gas and continuum 
dust emission. Our SMA observations resolve the disc structure 
down to scales of 200 AU, and allow us to probe the structure of 
the inner disc (< 400 AU) and the outer disc (400-900 AU). Our 
main conclusions can be summarized as follows. 

- The '^CO and '^CO emission extends to 900 AU from 
IM Lup, much further than the outer radius of 400 AU in- 
ferred earlier from dust measurements. 

- The H2 gas surface density in the region between 400 and 
900 AU Hes in the range of 5 x 10^" to 10^^ cm"^, using the 
standard CO-to-H2 ratio of 10 

- The disc is in Keplerian rotation around a central mass of 
1.2 + 0.4 Mq. Infall motions, if present in the outer disc, are 
negligible at < 0.2 km s"'. 

- The molecular line emission from the inner disc, withi n 
400 AU, is well described by the model of lPinte et all d2008l) . 
except that the gas temperature in the layers dominating the 
line emission of '^CO and '^CO exceeds the dust tempera- 
ture by factors 1 .7 and 1 .4, respectively. 

- Outside 400 AU, the surface densities of the molecular gas, 
as traced through '^CO and '^CO, of small (~ 1 pm) dust 
grains, and of larger (~ 1 mm) dust grains have a break in 
their radial dependence. At 400 AU, the dust surface density 
(in small grains) drops by a factor ~100, while the gas sur- 
face density shows a comparable drop of a factor 10-200 or 
steepens its radial power-law slope from p - lto5</7<8. 
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Our observations show that the disc around IM Lup consists 
of two regions. The inner 400 AU is well described by a 'stan- 
dard' accretion disc; the region between 400 and 900 AU has a 
much lower surface density as traced through dust grains with 
sizes from ~ 1 yum to 1 mm and through CO emission. Our ob- 
servations do not tell us if this outer region consists of material 
from the original prestellar core that has not (yet) made it onto 
the viscous accretion disc, or of material that is part of the disc 
but has had a different evolution. Sensitive, spatially resolved 
observations at various (sub) millimetre wavelengths, as may be 
obtained with the Atacama Large Millimeter Array may help to 
assess whether significantly different grain populations exist in- 
side and outside of 400 AU. With the same telescope, very high 
signal-to-noise observations of '-CO lines at high spectral res- 
olution may allow determination of any radial (inward or out- 
ward) motions in the > 400 AU gas. Spatially resolved mid- 
infrared imaging in several emission bands of PAHs, as could be 
obtained with the VISIR instrument on VLT, would shed light on 
the question if the 400-900 AU zone in the disc is largely pho- 
todissociated or -ionized. Detailed modelling of dust evolution 
in the outer disc may answer whether radial drift is responsible 
for the low column of dust beyond 400 AU in IM Lup disc. 
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